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We investigate the chiral and the deconfinement transition within the 
framework of Dyson- Schwinger equations using quenched lattice data for 
the temperature dependent gluon propagator as input. We extract corre- 
sponding order parameters from the Landau gauge quark propagator with 
?7(l)-valued boundary conditions. We study the chiral transition using the 
conventional quark condensate, whereas for the deconfinement transition 
we determine the dual condensate ('dressed Polyakov loop'). In addition we 
consider an alternative order parameter for deconfinement, the dual scalar 
quark dressing function. As a result we find almost the same transition 
temperatures for the chiral and deconfinement transitions. 

PACS numbers: 12.38.Aw, 12.38.Lg,11.10.Wx 

1. Introduction 

The phases of QCD are currently under intense theoretical and experi- 
mental investigations. Open questions concern among others the interplay 
between the confinement /deconfinement transition and the chiral transition, 
e.g. the (non-) coincidence of the chiral and the deconfinement transition at 
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zero chemical potential [I], and the possibility of a confined chirally sym- 
metric ('quarkyonic') phase [2]. Answers to these questions certainly require 
non-perturbative approaches to QCD. 

Functional methods involving the renormalization group equations [3j 
and/or Dyson-Schwinger equations (DSE) [4, 5j constitute a non-perturbative 
continuum approach to QCD. Only very recently, approaches accounting 
also for the deconfinement transition became available for these methods 
[BJ [U [9] . These approaches allow to define order parameters for the de- 
confinement transition from the quark propagator for generalized boundary 
conditions. Originally introduced within the lattice framework [10\ [TT] it 
was adapted to functional methods in El [9] . The quantities calculated 
in this work signaling the deconfinement transition are the dual quark con- 
densate (or 'dressed Polyakov loop') and the dual scalar quark dressing 
function. 

In the following we first recall the defining equations for the ordinary 
and the dual quark condensate and scalar dressing function, then summarize 
the truncation scheme used in our DSE calculations before we discuss our 
results for the chiral and deconfinement phase transition. 

2. Order parameters for deconfinement 

Consider the full quark propagator at finite temperature given in terms 
of its Dirac structure by 

S(p, lu p ) = [174 u; p C(p, uj p ) + i-fi PiA(p, u p ) + B(p, , (1) 

with vector dressing functions A and C and scalar dressing function B. 
The physical, antiperiodic boundary conditions in temporal direction lead 
to u>p(nt) = (27rT)(nj + 1/2) for Matsubara frequencies. Generalizing to 
temporal [/(l)-valued boundary conditions ip(x,l/T) = e ltp il)(x,ti) results 
in Matsubara frequencies u p (nt) = (2-7rT)(rij + (p/2) with boundary angle 
99=[0,2vr[. 

From the (^-dependent propagator the dual quark condensate introduced 
in lattice gauge theory Ref. [11] 

21 = 1 ^e-*6fy>„ (2) 
o 

can be calculated. Due to its close connection to the Polyakov loop it is 
also called 'dressed Polyakov loop', see Refs. [HI [12] for more details. 

The dual scalar quark dressing function introduced in Ref. [9] is the 
phase-Fourier-transform of the (/^-dependent scalar quark dressing function 
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Fig. 1. The Dyson- Schwinger equation for the quark propagator. Filled circles 
denote dressed propagators whereas the empty circle stands for the dressed quark- 
gluon vertex. 



evaluated at lowest Matsubara frequency and zero momentum 

Sb = | £(0,^(0,^)). (3) 

o 

Both quantities the dual quark condensate and the dual scalar quark dress- 
ing function transform under center symmetry identically as the conven- 
tional Polyakov loop and are therefore order parameters for deconfinement. 

An advantage of the 92-dependent quark condensate (ipip)^ is its direct 
connection to the ordinary Polyakov loop in the limit of static quarks 
For vanishing bare quark mass it is well-behaved in the continuum limit 
whereas at finite quark mass it is quadratically divergent and needs to be 
properly regularized. Though apparently not linked to the Polyakov loop the 
99-dependent scalar quark dressing function B(0, uj p (0, ip)) has the advantage 
of being a well defined quantity also in the continuum limit. 

3. The Dyson-Schwinger equation for the quark propagator at 

finite temperature 

The Dyson-Schwinger equation for the quark propagator is displayed 
diagrammatically in Fig. [TJ At finite temperature T it is given by 

S-\p) = Z 2 S^(p) - C F Z XJ g 2 T^ J —^S(k)T u (k,p)D flI/ (p-k)(4) 

with p = (p,uj p ) and k = (kjUJk) and renormalization factors Z2 and Z\f. 
Here denotes the (transverse) gluon propagator in Landau gauge and 
T u the quark-gluon vertex. The bare quark propagator is given by Sq 1 (p) = 
i'y-p+m. The Casimir factor Cf = (N^ — 1)/N c stems from the color trace; 
here we only consider the gauge group SU(2). The quark dressing functions 
A,B,C can be extracted from Eq. (JH) by suitable projections in Dirac-space. 

In order to solve this equation we have to specify explicit expressions for 
the gluon propagator and the quark-gluon vertex. At finite temperatures the 
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Fig. 2. Quenched SU(2) lattice results [T3] for the transverse dressing function 
Zx(q) and the longitudinal dressing function Z^(q) of the gluon propagator to- 
gether with the fit functions [S]. 



tensor structure of the gluon propagator contains two parts, one transversal 
and one longitudinal to the heat bath. The propagator is then given by 
(q = (q,UJ q )) 

D,M) = ^-P T ,M + ^ p iM (5) 

with transverse and longitudinal projectors 

p»M = (s^ - ^) Mi* , P^(q) = iVG?) - p^(q) , (6) 

with (i,j = 1...3). At zero temperatures Euclidean 0(4)-invariance re- 
quires both dressing functions to agree, i.e. Zx(q) = = Z{q). 

The temperature dependence of the gluon propagator can be inferred 
from recent lattice calculations. The results of Ref. [13] are shown in Fig. [2j 
The lattice data although still with quantitative uncertanties may very well 
correctly represent the qualitative temperature dependence of the gluon 
propagator. We therefore use a temperature dependent (qualitative) fit to 
the data as input into the DSE; this fit is also displayed in Fig. [2] (straight 
lines). The fit functions are described in detail in Refs. [3(9] and shall not 
be repeated here for brevity. Note, however, that we also inherit the scale 
determined on the lattice using the string tension y/a = 0.44 GeV |13| . 

For the quark-gluon vertex we employ a temperature dependent model 
which is discussed in detail in [9]. 
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Fig. 3. Left diagram: Temperature dependence of the dressed Polyakov loop Si 
and the conventional quark condensate A„ = {ip^) v= Tt together with their deriva- 
tives for m = lOMeV. Right diagram: The same quantities in the chiral limit. 



4. Numerical results 

In Fig. [3] we display our numerical results for the ordinary and the dual 
quark condensate together with their (normalized) temperature derivatives 
once evaluated for a quark mass of m = lOMeV (fixed at T = 200 MeV 
and /3 2 = 20GeV 2 ) and once evaluated in the chiral limit. One clearly sees 
the difference in the chiral transition: whereas at finite bare quark mass 
we encounter a crossover the transition changes into a second order phase 
transition in the chiral limit. In the first case the corresponding tempera- 
ture derivative shows a peak at T c = 301(2) MeV, whereas it diverges at 
T c = 298(1) MeV in the second case. We also extracted the corresponding 
transition temperatures from the chiral susceptibility 

XR = m2 -^ (W)t - (V^)t=o) • ( 7 ) 

The results for quark mass m = 10 MeV are given in table [TJ 

The corresponding transition temperature for the deconfinement transi- 
tion can be read off the dual quark condensate (or dressed Polyakov loop). 
At finite quark mass and in the chiral limit we observe a distinct rise in 
the dual condensate around T ~ 300 MeV. The corresponding (normalized) 











301(2) 


304(1) 


305(1) 


308(2) 



Table 1. Transition temperatures for the chiral and deconfinement transition for 
quark mass to = 10 MeV. 
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Fig. 4. Left diagram: Comparison of the temperature dependence of the corre- 
sponding dual scalar quark dressing and the dual condensate for m — 10 MeV. 
Right diagram: Temperature dependence of the chiral quark condensate A v = 
{ipi)}ip and the scalar dressing function B v at periodic boundary conditions ip = 0. 



temperature derivative shows a peak at Td ec = 308(2) MeV for quark mass 
m = 10 MeV. In the chiral limit this peak moves to T^ ec = 299(3) MeV. 

In general we note that the chiral and deconfinement transition are close 
together. There are a few MeV difference between the different transition 
temperatures for the crossover at finite quark masses, whereas both tran- 
sitions occur at the same temperature (within errors) in the chiral limit. 
These findings agree with early expectations from lattice simulations [14J. 

Furthermore we wish to emphasize that the present calculation, although 
carried out with quenched lattice results for the gluon propagator, is in 
itself not strictly quenched. This can be seen from the fact that the dressed 
Polyakov loop is not strictly zero below the deconfinement transition. We 
refer here to quenched calculation since the influence of the matter content 
to the gauge sector is neglected. 

In the left diagram of Fig. 0] we show a comparison of the dual scalar 
dressing Eq. ([3]) with the dual condensate for m = 10 MeV. Both quantities 
show a very similar temperature dependence especially the aforementioned 
distinct rise around T ~ 300 MeV. The signal is slightly stronger for the 
dual scalar dressing S b ■ The temperature derivatives indicating the decon- 
finement transition clearly peak at the same transition temperature. 

In the plot of the right hand side of Fig. 0] the chiral condensate (V^^o 
and the quark scalar dressing B(ip = 0) as functions of temperature at pe- 
riodic boundary conditions are shown. For comparison both are normalized 
to 1 at T = 500 MeV. Whereas (V>V) v=0 is a strictly monotonic function 
with temperature the scalar dressing is slightly decreasing in the tempera- 
ture range between 300 and 400 MeV before the high temperature behavior 
shows up. The large temperature scaling of both quantities can be extracted 
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analytically from Eqs. ([2]) and (Jl|) as pointed out in the appendix of Ref. [9]. 
There it is shown that 

^ =0 , P=0 (T) ~ VT for T > T c (8) 

and consequently a quadratic rise of the condensate, 

(^=0 ~ T 2 , for T » T c , (9) 

is obtained. These results are in excellent agreement with polynomial fits to 
the data. The dual condensate Si as well as the dual scalar dressing £g are 
promising candidates for a further study of the deconfinement transition of 
QCD. 

5. Summary 

In this talk we addressed the chiral and the deconfinement transition 
of quenched QCD. We showed results for the order parameter of the chiral 
transition, the quark condensate, and for order parameters of the decon- 
finement transition, the dressed Polyakov loop and the dual scalar dressing. 
These were extracted from the Landau gauge quark propagator evaluated at 
a continuous range of boundary conditions for the quark fields. Independent 
of the dual order parameter used we found exactly the same deconfinement 
transition temperature. A comparison of the transition temperatures for 
the chiral and the deconfinement transition shows almost coincidence for 
moderate quark mass of the order of an up-quark. In the chiral limit the 
two transitions coincide within error. We find a second order chiral phase 
transition at T Xr j t a = 298(1) MeV and a similar temperature for the decon- 
finement transition, Td ec = 299(3) MeV. In summary we conclude that the 
chiral and deconfinement transition temperatures are only slightly different 
for finite quark masses and coincide within errors in the chiral limit. 

The framework used in this work is quenched SU(2) Yang-Mills the- 
ory. Our transition temperature may be translated into the corresponding 
ones of quenched 5*7(3) QCD using the relations T c /^/a = 0.709 (577(2)) 
and T c l yfa = 0.646 (577(3)) between the respective critical temperatures 
and the string tension [IB]. The resulting transition temperature is then 
T XR / T 4 Rj Tdec ~ 272 MeV in the chiral limit. In order to work in the 
full, unquenched theory quark loop effects and meson effects which shift 
the transition temperature below T = 200 MeV have to be taken into ac- 
count. Concerning the dual condensate and scalar dressing function in the 
unquenched formulation additional effects due to the Roberge- Weiss sym- 
metry |15] occur. This is because of the formal similarity of the continuous 
boundary conditions for the quark field to an imaginary chemical potential, 
see [6J for details. 
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